Numerical study of hydrogenic effective mass theory for an 
impurity P donor in Si in the presence of an electric field and interfaces 
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In this paper we examine the effects of varying several experimental parameters in the Kane 
quantum computer architecture: A-gate voltage, the qubit depth below the silicon oxide barrier, 
and the back gate depth to explore how these variables affect the electron density of the donor 
electron. In particular, we calculate the resonance frequency of the donor nuclei as a function of 
these parameters. To do this we calculated the donor electron wave function variationally using an 
effective mass Hamiltonian approach, using a basis of deformed hydrogenic orbitals. This approach 
was then extended to include the electric field Hamiltonian and the silicon host geometry. We 
found that the phosphorous donor electron wave function was very sensitive to all the experimental 
variables studied in our work, and thus to optimise the operation of these devices it is necessary to 
control all parameters varied in this paper. 

PACS numbers: 03.67.Lx, 71.55.Cn, 85.30.De 
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I. 



INTRODUCTION 



Since Kohn and Luttinger'si*^ original work on shal- 
low donors in silicon, there has been renewed interest 
in the study of donor impurities in silicon, particularly 
the Si:"^^P system, following Kane's^ proposal for a solid- 
state quantum computer. In the Kane quantum com- 
puter, information is encoded onto the nuclear spins of 
donor phosphorous atoms in doped silicon electronic de- 
vices. Application of an electrostatic potential at surface 
electrodes positioned above the qubits (A-Gates) tunes 
the resonance frequency of individual spins, while surface 
electrodes between qubits (J-Gates) induces electron- 
mediated coupling between nuclear spins. Perturbing the 
donor electron density with an externally applied elec- 
tric field is crucial in tuning the hyperfine interaction 
between the donor electron and nucleus and hence also 
in tuning the resonance frequency of the P nuclei and 
controlling logical operations. Substantial theoretical ef- 
forts have been devoted to modeling the P donor electron 
ground state in the silicon wafer device, and the altered 
ground state with an externally applied electric field. In 
this paper we discuss relevant experimental parameters 
which can be controlled to perturb the donor electron 
wave function. 

There is a considerable amount of work done in this 
area, and several theoretical approaches have been pur- 
sued with varying degrees of application and approxi- 
mation. In Kohn and Luttinger's^'^ work, the P donor 
ground state in the bulk silicon is calculated using a sin- 
gle trial wave function: a deformed IS* hydrogenic orbital 
and varying the Bohr radii to minimise the ground state 
energy. In this paper we follow Faulkner's^ approach and 
extend Kohn's method to include a trial wave function 
expanded in a basis of deformed hydrogenic orbitals, and 



vary the Bohr radii to minimise the ground state energy. 
As we have used a large basis set in this approach, the 
ground state wave function has the flexibility to distort 
with the application of an electric field above the P donor. 
Several authors^iSii have previously investigated the ef- 
fects induced by strain and interface regions on donor 
states. These external influences partially lift the valley 
degeneracy in the bulk silicon. 

The effect of an electric field potential at a gate above a 
P donor in a silicon substrate on the hyperfine interaction 
coupling between the P donor electron and nucleus has 
already been reported by several authors. In the work 
of Kane^ and Larinov et al.,^ the effect of an electric 
field potential in the bulk silicon host is considered using 
perturbative theory, excluding the additional interface 
potentials. Wellard et alM> consider both the influence 
of the electric field and interface barriers using a spherical 
effective mass Hamiltonian. 

The main advantage demonstrated in our approach us- 
ing the anisotropic basis is the flexibility in choosing the 
smaller effective Bohr radius for the donor ground state 
to be in the direction towards the interface regions. This 
minimizes the overlap of the donor wave function into 
these regions. For shallow donor depths, the donor wave 
function is restricted in moving towards the A-gaXe be- 
cause of the silicon oxide interface. 

In this work, we include the effects of both the elec- 
tric fleld potential and the interface regions, and the 
anisotropy of the conduction band minimum in Si. To 
our knowledge, there have been hitherto no published 
results for modeling electrostatic gate operations in the 
Kane quantum computer which include simultaneously 
the anisotropy of the effective masses in the silicon host, 
the electric fleld potential and the interface regions in the 
Si wafer device. In this paper we address all these criteria 
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and discuss relevant experimental parameters which can 
be adjusted to fine tune the contact hyperfine interaction. 
We calculate this coupling as a function of ^-gate volt- 
age, donor depth below associated ^-gate and the back 
gate depth. A subsequent paper will discuss our further 
results for the J-gate controlled electron exchange inter- 
action between adjacent donor electrons. 

In Sec. ini we will discuss some background effective 
mass theory and the approach we took to obtain the 
phosphorous donor ground state in bulk silicon with no 
electric field applied. Section IIIII discusses how we ob- 
tained the electric field potential and modeled the sili- 
con host geometry to include the silicon oxide layer and 
back gate. The numerical results using the methods out- 
lined in the previous sections are presented in Sec. II VI for 
the varying experimental parameters studied. Finally we 
summarise our major findings in Sec. IVl 
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FIG. 1: Ground state electron density without electric field. 



II. 



FAULKNER'S METHOD 



19.94meV, Eq. becomes: 



Neglecting inter-valley terms, the one-valley effective 
mass equation for the energy levels of donors in silicon is 
given below 1^ 



2m± \dx'^ dy"^ 



92 



2m II dz^ 



e 

er 



5f(r) = E^{r), 
(1) 



where e = 11.4 is the dielectric constant, and m± = 
0.1905mo and m|| = 0.9163mo are the transverse and 
longitudinal effective masses respectively, and toq is the 
mass of a free electron. Here we are expanding the energy, 
E'^ around the conduction band minimum along the z- 
axis at A; = (0, 0, fco)- 



2ra 



2m\\ 



kl). (2) 



We followed Faulkner's approach and kept the full 
anisotropy of the conduction band minimum. We ex- 
panded the donor electron wave function, ^'(7*) in a basis 
of deformed hydrogenic orbitals: 



1/4 



■>pnim{x,y, J —z,a), 



where 'ipni^{x,y, z, a) = Rni{a,r)Yira{0,4>) are the nor- 
malised hydrogenic orbitals, 7 — m±/m\\ ~ 0.2079, a is 
the effective Bohr radius in the x,y directions, and (3 is 
an adjustable parameter which gives the effective Bohr 
radius h in the z direction. 

If we use atomic units, where the unit of length 0^ — 
S^e/TOj_e^ = 31. 7A and unit of energy 'm^e^/2'h?e^ = 



92 
dy- 



+ 7 



dz^ 



*(r) = £;*(r) 



(3) 



Equation Q was diagonalised with the effective Bohr 
radius a and parameter (3 varied to minimise the ground 
state energy E. The ground state energy converged using 
a basis of 91 hydrogenic orbitals to give E — — 3 1.23m eV, 
and effective Bohr radii: a = 23.8lA and b = \/ (7//3)a — 
13.68A. These results are consistent with Kohn's results 
of a = 25A and b = 14.2,^1*^ and Faulkner's ground 
state energy E — — 31.27me\4i for phosphorous. 

The ground state wave function obtained was a de- 
formed hydrogenic IS' orbital. Figurc^shows the ground 
state electron density plotted in the x, y and z direc- 
tions for comparison of the different effective Bohr radii 
obtained in the different directions. Also shown in this 
figure is the ground state obtained using a spherical effec- 
tive mass Hamiltonian and isotropic hydrogenic orbitals 
as a basis, here the effective electron mass is given by 
TO* w TO_L = 0.1905mo, which gives an effective Bohr 
radius of 3.17nm. 



III. INCLUDING THE ELECTRIC FIELD AND 
SILICON HOST POTENTIAL 

Faulkner's method was then extended to include the 
effects of an electric field above the qubit, and boundary 
conditions of the silicon host. The solution of Poisson's 
equation to extract the electric field potential for our de- 
vice with the A-gate at varying voltages was obtained by 
simulation using a Technology Computer Aided Design 
(TCAD) modeling package. 

TCAD is used in the electronics industry as a tool for 
2-D and 3-D modeling and simulation of semiconductor 
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FIG. 2: Schematic design parameters implemented in TCAD 
to model the Kane computer architecture. 



devices. It employs a coupled Newton-like solver at dis- 
crete nodes to obtain the self-consistent solution of the 
Poisson and electron- hole continuity equations. Figure El 
shows the 2-D device scheme implemented in TCAD used 
to model the application of a voltage to the A-gate above 
qubit, Qi. The lateral edges of the silicon lattice were 
assumed to extend infinitely in the y-direction, but the 
electrostatic potential was only obtained on a finite grid 
210nm wide, with the potential set to zero outside this 
region. We checked that this approximation is valid at 
the boundaries and found the TCAD potential had fallen 
close to zero (10~^ — 10~^eV), at ?/ = ±105nm. The po- 
tential in 2-D from TCAD is assumed to have a "thick- 
ness" in the third dimension (x) of Ifim. 

In this paper we examine the effects of varying several 
experimental parameters: A-gate voltage, qubit depth 
below the silicon oxide barrier, and the back gate depth 
to explore how these variables affect the electron density 
of the donor electron at the phosphorous nuclei. In par- 
ticular we calculate the resonance frequency of the donor 
nucleus as a function of these parameters. 

The application of a potential, and the silicon host ge- 
ometry in the device shown in Fig.[5]splits the degeneracy 
of the two local minima along the z-axis, compared to the 
other four along the x and y axis in the lower Ai ground 
state. ^ With no electric field applied the ground state 
wave function is > 99% 15* in character. When the volt- 
age applied is low enough so that the wavefunction stays 
predominantly IS' in character, diagonalising the single 
valley effective mass equation is equivalent for solving in 
either valley, ±z, since the deformed IS* wavefunction is 
symmetric in z«S. 

Using these justifications we can formulate the problem 
using a co-ordinate system with the z-axis in the direc- 
tion from Qi to the interface. With this convention we 
expand the donor wave function around the conduction 
band minimum oriented along the z-axis. Because of the 
smaller effective Bohr radius in the z direction towards 



the interface and back gate, the ground state is lower in 
energy since there is less penetration of the wave function 
into these barrier regions. 

With the electric field the Hamiltonian is: H = 
Hq + iJi, where Hq is the zero field Hamiltonian, and 
Hi ~ V{y, z) is the electric field potential term. V{y, z) 
is the electric field potential generated from TCAD, and 
here we also add an additional term to model the Si02 
layer and the back gate. The Si/Si02 barrier was mod- 
eled as a step function with height 3.25eV, since most 
insulators have a work function greater than 3eV>ii The 
back gate serves as a reference voltage point (ground) to 
the voltages applied to the top gates. Outside the back 
gate the potential was set at 3.25eV also. 

To calculate the perturbed donor electron wave func- 
tion and energies we constructed the electric field Hamil- 
tonian matrix. Hi , with its elements given by: 



{n'l'm'\Hi\nlm) 



y-i^nimix, y, \l -z, a). 



(4) 



The integrals in Eq. @ were then calculated numer- 
ically for the varying voltages at the A-gate and qubit 
position. Once Hi was obtained the total Hamiltonian 
was then diagonalised to find the donor electron ground 
state with the varying experimental parameters. 



IV. 



NUMERICAL RESULTS 



A. Results obtained varying A-gate voltage and 
donor depth 

The perturbed donor electron ground state was cal- 
culated for each set of experimental parameters to com- 
pare and optimise the conditions for addressing the tar- 
get qubit, Qi. Once the electron ground state was found 
we calculated the value of the contact hyperfine coupling 
A{V) for each particular voltage at the A-gate and qubit 
depth below this gate. 

The general formula for the contact hyperfine coupling 
A{V) is given below: 

MV) = ^A^sWw/iol*("t^,0)p, (5) 

where ^(V,0) is the donor electron ground state wave 
function evaluated at the donor nucleus, /x_b is the Bohr 
magneton, is Lande's factor for ^^P, /xjv is the nuclear 
magneton and /xq is the permeability of free space»2ii£ 

Since we use effective mass theory, instead of calculat- 
ing the donor wave function with the full expansion of 
the Bloch functions, we calculate the envelope function. 



4 




following equation (accurate to second order)* 
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FIG. 3: Nuclear resonant frequency shifts of qubit Qi at 
lower voltages with varying donor depths, back gate depth at 
60nm, using anisotropic basis. 
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FIG. 4; Nuclear resonant frequency shifts of qubit Qi at 
varying voltage and donor depth, back gate depth at 60nm, 
using anisotropic basis. 



which describes the smooth donor-related modulation of 
the electron wave-function. So instead of calculating the 
contact hyperfine coupling, A{V), directly we calculate 
the relative shift in A{V) with the potential applied and 
assume this shift will be similar to those of the true wave 
functionJiS Thus we need to calculate: 



AiV) 



(6) 



where A{0)/h = 28.76MHz is determined for ^^P in sili- 
con from experimental dataj^ and ^'(V, r) are the donor 
envelope wave functions calculated by our method. 

The phosphorous nuclear resonant frequency is affected 
by the donor electron when the valence electron is spin 
polarized by a background magnetic field, B, of the order 
of 2T. The hyperfine interaction constant is related to 
the frequency separation of the nuclear levels, via the 



hv 



2g]ypL]yB + 2A- 



2A^ 



(7) 



In all the calculations we considered the background 
magnetic field fixed at 2T. Figure 01 shows the nuclear 
resonant frequency shift of Qi, calculated for a lower 
range of positive A-gatc voltages, between OV and 0.8V, 
for the varying donor depths below the silicon oxide bar- 
rier. Figure^ shows the nuclear resonant frequency shift 
calculated for the full range of j4-gate voltages, between 
— l.OV and l.OV, for the varying donor depths below the 
silicon oxide barrier. These plots are calculated with a 
close back gate depth set at 60nm 

For comparison of our method with previous resultsifi 
reported using a spherical effective mass Hamiltonian, 
we calculated the resonance frequency of Qi using an 
isotropic effective Bohr radius of « 3nm. Our results 
were consistent with the calculations of Wellard et al)^ 
The results for the isotropic basis showed that for donor 
depths close to the silicon oxide barrier, the wave func- 
tion was restricted in moving towards the applied A-gate 
voltage. The donor wave function obtained using the 
anisotropic basis, is advantageous because of the smaller 
effective Bohr radius in the direction toward the sili- 
con oxide layer, which results in less penetration of the 
donor wave function into the interface regions. Thus the 
anisotropic basis produced a more energetically favorable 
ground state than the isotropic ground state. 

For the lower voltages (< 0.8V), the results are con- 
sistent with the expectation that the closer the donor 
depths are to the applied voltage, the greater the fre- 
quency shift. At voltages above a certain threshold and 
donor depths further away from the silicon oxide bar- 
rier, there is a huge difference in the donor wave function 
from the zero field ground state, as it is perturbed almost 
completely away from the nucleus. Figure shows an ex- 
ample of this change in electron density for a voltage of 
l.OV at the A-g&tc and donor depth of 40nm. Here the 
P nucleus is at the origin and as z decreases the electric 
field increases. 

In Fig. 13 and [7| we observe the difference in the donor 
electron ground state obtained for a donor depth of 20 
and 40nm with a positive voltage of l.OV at the A-gate. 
In both these plots the donor wave function moves to- 
ward the applied A-gate voltage in the negative z direc- 
tion. For a close donor depth of 20nm we observe that 
even though the donor wave function moves slightly to- 
ward the A-gate, it is significantly restricted in moving in 
this direction because of the silicon oxide interface in this 
direction also. In contrast the donor wave function for 
a depth of 40nm deforms unhindered toward the A-gate, 
and most of the electron density has been transformed 
away from the nucleus. 

Figure |S| shows the donor electron density obtained in 
the yz-plane for a negative voltage of -l.OV at the A- 
gate and a donor depth of 20nm. A negative applied 
voltage causes the electron to disperse in all directions 
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FIG. 5: Ground state electron density in yz-plane for donor 
depth at 40nm and voltage at IV at the A-gate. 
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FIG. 6: Ground state electron density for donor depth at 
20nm and voltage at 0. and l.OV at the A-gate, in the z- 
direction. 
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FIG. 7: Ground state electron density for donor depth at 
40nm and voltage at l.OV at the A-gate, in the z-direction. 
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FIG. 8: Ground state electron density in yz-plane for donor 
depth at 20nm and voltage at -l.OV at the yl-gate. 

away from the positive potential, this plot demonstrates 
that because of the close back gate in the positive z- 
direction, the electron density predominantly perturbs 
away from the applied voltage in either direction laterally. 

Because of the interface regions it is either energeti- 
cally favorable for the donor electron wave function at 
shallow donor depths to distort completely away from 
the nucleus, when the gate voltage is negative, or for the 
donor wave function to be restricted in distorting towards 
the ^-gate, with a positive voltage. 

In Table HI we present a comparison of the difference 
in the ground state energy for the donor wave function 
without the electric field (Eq), and with a positive voltage 
of l.OV apphed to the A-gate Also reported in 

this table is the TCAD potential at the P nucleus for the 
varying donor depths. 



TABLE I: Eiv — Eq for a back gate depth of 60nni. 



Qi Depth 


TCAD Potential 
at Qi (meV) 


Eiv — Eo 
(meV) 


20nm 


-90.02 


-91.70 


40nm 


-37.06 


-47.73 



For the close donor depth at 20nm we observe that the 
energy difference is approximately equal in magnitude 
to the TCAD potential at the nucleus. This is because 
the donor wave function has perturbed only slightly from 
the zero field ground state wave function. In contrast the 
energy difference for the donor depth at 40nm is much 
higher as the wave function deforms significantly from the 
ground state wave function towards the applied voltage. 

If we compare the results obtained in our work includ- 
ing the effect of the interface barriers in addition to the 
electric field potential, with Kane's® results wherein only 
the potential of a uniform electric field in the bulk was 
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considered, we observe that the siHcon oxide layer and 
the back gate exert a significant influence on the donor 
electron's ground state. Instead of the contact hyperfine 
coupling, A{V), being independent of whether a positive 
or negative voltage is being applied at the A-gate as re- 
ported by Kane, we observe in Fig. 01 that the interface 
regions in the silicon host geometry break this symmetry. 

Even without considering the influence of the interface 
regions, the effect of whether a positive or negative volt- 
age is applied at the A-gate causes very different changes 
in the donor electron density. For a positive voltage the 
electron is bound to both the nucleus and the A-gate. In 
contrast, when a high enough negative voltage is applied 
so that the electron is no longer bound to the P nucleus, 
the electric field profile causes the electron to disperse in 
all directions away from the positive potential. 



B. Results obtained varying back gate depth and 
donor depth 
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FIG. 9: Nuclear resonant frequency shifts of qubit Qi at 
varying donor depths with back gate depth at 60 and lOOnm, 
and l.OV at A-gate. 



To observe the effect that the back gate depth has on 
the donor electron wave function we repeated the calcu- 
lation with a back gate depth at lOOnm. Figure |H| shows 
the comparison between nuclear resonant frequency shifts 
of the donor electron with the application of a voltage at 
the A-gate with a close and far back gate. These calcu- 
lations were performed with a close back gate at 60nm 
and a far back gate at lOOnm, with a bias of l.OV at the 
A-gate and donor depths ranging from 30 to 75nm. 

With a closer back gate the electric field strength was 
higher within the Si wafer, and the donor electron wave 
function was perturbed greater, and so the frequency 
shift was more pronounced for donor depths with a close 
back gate. For donor depths close to the back gate the 
interface barrier effectively "pushes" the electron towards 
the A-gate. With the back gate at lOOnm, the electric 
field strength is lower, and there is no substantial overlap 
of the donor electron wave function with the back gate 
barrier for donor depths of 30 and 40nm, so it is not as 
energetically favorable for the donor electron to perturb 
away from the back gate toward the A-gate. 

Figure^lshows the ground state wave function plotted 
in the yz-plane for a donor depth of 75nm and with a 
back gate depth of lOOnm and a positive voltage of l.OV 
at the A-gate. This plot demonstrates that even at a 
donor depth far from the A-gate, the ground state wave 
function distorts freely toward the A-gate because of the 
close proximity of the back gate, and the remoteness of 
the silicon oxide interface. 

In Table HTl we present a comparison of the difference in 
the ground state energy for the donor wave function with- 
out the electric field, and with a voltage of l.OV applied 
to the A-gate, and the back gate at lOOnm. Also reported 
in this table is the TCAD potential at the P nucleus for 
the varying donor depths. This table reflects the trend 
noted in Table ^ that a signiflcantly lower ground state 
energy is obtained for the deeper donor depths, where 



m- X lo'^ 

5 r 




FIG. 10: Ground state electron density in j/z-plane for donor 
depth at 75nm with back gate depth of lOOnm, and l.OV at 
A-gate. 



the electron density perturbs significantly away from the 
nucleus toward the applied voltage. 



TABLE II: Eiv - Ea for a back gate depth of lOOnm. 



Qi Depth 


TCAD Potential 
at Qi (meV) 


Eiv — Eo 
(meV) 


40nm 


-67.25 


-68.33 


75nm 


-26.00 


-36.57 
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V. CONCLUSIONS AND IMPLICATIONS FOR 
CURRENT FABRICATION TECHNOLOGY AND 
DEVICE MODELING 

We believe that the resuhs reported here using effective 
mass theory are quantitatively reasonable. It is reason- 
able to expect that the variation of the donor wave func- 
tion with the experimental parameters calculated here 
using the smooth donor envelope function, would be sim- 
ilar to that of the true wave function. 

It is evident that the P donor electron wave function is 
sensitive to all experimental parameters studied in this 
paper. The donor wave function exhibits a fundamen- 
tal change at crucial experimental parameters, where 
the electron wave function transforms from being only 
slightly perturbed from the zero field ground state, to be- 
ing almost c;omplc;tcly perturbed from the nucleus. These 
results highlight the significance of the influence of the 
silicon host geometry on the donor electron wave func- 
tion. Ongoing work in our laboratory is focusing on ver- 
ification of these results, using the full Bloch wave struc- 
ture in our calculations. These results demonstrate the 
importance of the boimdary conditions imposed by the 
interface regions, and the need to use a basis set which 



has the flexibility to meet the boundary conditions. 

However, including the Bloch wave structure, the inter 
valley terms and the electric field and interface potentials 
is a challenging task. The results presented are quantita- 
tively reasonable and provide a fast and reliable method 
which gives insight into the behavior of the P donor elec- 
tron wave function under several different experimental 
conditions. To optimise the fabrication of these devices 
it is necessary to take into account the dependence of 
the donor electron wave function on all parameters var- 
ied in this paper: donor depth below the A-g&te, back 
gate depth and voltage at the ^-gate. 
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